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We consider the semiclassical approximation to the spectral form factor K{t) for two- 
dimensional uniformly hyperbolic systems, and derive the first off-diagonal correction for small 
r. The result agrees with the r^-term of the form factor for the GOE random matrix ensemble. 
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A fundamental characteristic of quantum systems that are chaotic in their classical limit 
is universality. It is observed that diverse systems behave identically when statistics of energy 
levels or wave functions are considered, provided that they have the same symmetries. These 
universal statistics agree with those of random matrix theory, i.e. with the statistics of eigenval- 
ues and eigenvectors of large random matrices |1| . Support for this random matrix hypothesis 
comes from a large number of numerical and experimental investigations which have been car- 
ried out on a great variety of systems . However, it remains an open question to understand 
the origin of this universality, and its relation to the underlying classical dynamics. 

One theoretical approach by which such an understanding may be attempted is the semi- 
classical method. Semiclassical approximations are asymptotically valid in the limit /?. — > 
where universality is expected to hold. Moreover, they directly connect quantum properties 
with properties of the corresponding chaotic classical system. They have been applied in par- 
ticular to statistical distributions of the energy levels which are bilinear in the density of states, 
one example being the spectral form factor K{t). One of the successes of the semiclassical 
approach has been to show that the spectral statistics do indeed agree with the random matrix 
statistics in the limit of long-range correlations; specifically the correct leading order behaviour 
of K{t) as r ^ has been derived 0. 

An extension of this result requires knowledge of correlations between different periodic 
orbits 0]. The relevant mechanisms by which periodic orbits are correlated have to be identified, 
and the contributions of correlated orbits to the spectral form factor have to be evaluated. 
Based on an analogy with disordered systems and with diffractive corrections it has been 
suggested that the next term in the expansion of K{t) for small r originates from 'two-loop 
orbits': orbits that have a self-intersection with small crossing angle and neighbouring orbits 
without self-intersection 0. There is strong numerical evidence that in systems with time- 
reversal symmetry these orbit pairs indeed yield the next-order-term in agreement with the 
expectation based on random matrix theory. 

In the following we present a derivation of the next to leading order term in the expansion of 
the spectral form factor for small r. We evaluate analytically the contributions of the two-loop 
orbits to the form factor for uniformly hyperbolic systems with time-reversal symmetry, and 
we show that the result indeed agrees with random matrix theory. The calculation makes clear 
the properties of classical trajectories which are responsible for the universal result. 

We consider the spectral form factor, which is defined as the Fourier transform of the two- 
point correlation function of the density of states 

where the density of states d{E) = J2n ^(-^ ~ ^n) is divided into a mean part d{E) and an 
oscillatory part do{E). For systems with time- reversal symmetry, or more generally an anti- 
unitary symmetry, it is expected that the form factor agrees in the semiclassical limit {h 0) 
with that of the Gaussian Orthogonal Ensemble (GOE) of random matrix theory which has 
the expansion 

K^o^{t) = 2t- 2t^ + 0(r3) as r ^ . (2) 

The semiclassical approximation for the form factor is obtained by inserting Gutzwiller's 
trace formula for the density of states into and evaluating the integral in leading order of h. 
The result is an approximation in terms of a double sum over all periodic orbits of the classical 
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Figure 1: The pairs of orbits considered here consist of different segments. In each segment one 
orbit is very close to the other (or its time reverse), but they can differ in the way the segments 
are connected. 

system 

*-m«s^(e^.^;-'-^^'<^-^)}^ (3) 

where r = T/{hd{E)) and h = 2TTh. Furthermore, is an amphtude, generally complex- 
valued, which depends on the stability and the Maslov index of the periodic orbit 7, and Sj 
and are its action and period. 

The double sum in (|^) runs over all possible pairings of periodic orbits. However, most of 
these pairs do not contribute in the semiclassical limit. Periodic orbits which are located in 
different regions in phase space are uncorrelated, and when summed over, the contributions from 
different pairs cancel each other. It is expected that the relevant semiclassical contributions 
come from a relatively small number of pairs of orbits which are correlated. The key problem 
is then to identify the mechanism which is behind these correlations. 

The basic assumption we make is that only those periodic orbits which are almost everywhere 
close to one another, or to the time-reverse of the other orbit, are correlated In order for 
two orbits to be different but nevertheless close they must have special forms which can be 
constructed in the following way. The orbits are composed of different segments during which 
one orbit follows very closely the other orbit (or its time-reverse) . However, the orbits can differ 
in the way in which the segments are connected. 

The two simplest possibilities are shown in Fig. 0. If orbits are composed of only one segment 
then the two ends can be connected in only one way. It then follows that the two neighbouring 
orbits are either identical or one is the time-reverse of the other. Including only these pairs in 
the double sum corresponds to the diagonal approximation, which yields the correct leading 
order behaviour K{t) ~ 2r as r — * 0]. 

Two segments, on the other hand, can be connected in two ways, leading to orbits with or 
without self- intersection at the connection point, as shown in Fig. |l|. In order for these pairs to 
exist and to be close, the crossing angle e has to be small. Then it can be shown in a linearized 
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approximation, that one orbit is indeed in the neighbourhood of the other. 

In the following we evaluate the contributions of such pairs of orbits to the spectral form 
factor. In order to avoid further assumptions and to keep the calculations simple we restrict 
attention now to systems with uniformly hyperbolic dynamics, specifically we consider the 
representative example of the geodesic motion on Riemann surfaces with constant negative 
curvature P]. Then the quantities A^, and Ty in (H) depend only on the length of an 
orbit, and Aj is positive. We assume that the systems have no further symmetries and are 
non-arithmetic so that the typical degeneracy of a length of a periodic orbit is two. For these 
systems the action difference for the pairs of orbits being considered here has been calculated 
in the linearized approximation for small crossing angle e in and is given by 



where A is the Lyapunov exponent of the system, and p and m are momentum and mass of the 
particle. 

The sum over these pairs of orbits can be evaluated by summing over all self-intersections 
of periodic orbits with small crossing angle e, because for every such self-intersection there 
exists a neighbouring periodic orbit with action difference AS{6). The self-intersections are 
determined by introducing a function which selects them. This is done in the following way. 
A self-intersection of a periodic orbit with period T divides the orbit into two loops. It can be 
characterized by the crossing angle e and the total time t along the shorter of the two loops, 
t < T/2. Furthermore, we introduce an angle variable that specifies the direction of the 
velocity, and a variable t' that measures the time along a periodic orbit. If at any time t' along 
a periodic orbit q(t' + t) = q{t') and (pit' + 1) = (pit') — ti + e then this periodic orbit has a 
self-intersection with opening angle e, and traversing the corresponding shorter loop takes time 
t. 

Correspondingly, we can express the contribution from pairs of the two-loop orbits to the 
form factor as 



delta functions to the three integration variables, where v is the speed of the particle and g is 
the determinant of the metric tensor. The three integrals give a contribution each time that t' 
is at the beginning of a loop with time t and opening angle e. The choice of the limits of the 
integral over e is not important since the main contribution in the semiclassical limit ^ — > 
comes from the asymptotic behaviour of the integrand at e = 0. In the amplitudes and the 
periods of the neighbouring orbits were set to be equal since the difference does not contribute 
to the leading semiclassical order. 

One of the important properties of long periodic orbits is their uniform distribution on the 
energy surface in phase space. It implies that the average of a given phase space function 




(4) 




(5) 



where the function f^.t is given by 
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f{q,p) along all periodic orbits of a certain period T can be replaced, in the limit T oo, 
by an average of this function over the energy surface in phase space More accurately, the 
following asymptotic relation holds as T ^ oo 



E,\A,mT-T,) J^dt' f{q{t'),p{t')) 

Jd'qd'p6[E-g^) f{q,p) (7) 



S{i?) 



where ^{E) is the volume of the energy surface in phase space. 

The relation (1) is in the form in which this property of the periodic orbits can be applied. 
The semiclassical limit ^ — > is performed with the condition that r/h ^ oo. The mean 
density of states being d{E) ~ T,{E)/{27ih)'^, this implies that T — > oo and thus the leading 
order semiclassical behaviour arises from the large T behaviour. Applying the uniformity of 
the periodic orbit distribution and performing the integral over the energy delta-function one 
obtains 

K^^\t) ^--^Re / dee^sm\e\ / dtpE{e,t) (8) 

mhd{E) J_„ Jo 

where 

PE{e, t) = j 5{q{t) - q,) - 0o + vr - (9) 

and q{t) and are the coordinates of a particle at time t, whose initial conditions at t = 
are specified by q^, (po and energy E. 

The quantity pE{£,t) has a direct classical interpretation. It is the probability density for a 
particle with energy E to return after time t to its starting point with a velocity that deviates 
from the initial velocity by an angle e — n. In the same way as for the diagonal approximation, 
one thus finds that the periodic orbit sum is related to a transition probability density in phase 
space 
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Our aim is to determine the leading order behaviour of (H) as h ^ which, as remarked 
above, depends on the long-time behaviour of pE{£,t). For long times pE{£,t) approaches one 
over the volume of the energy shell in phase space, because the particle is equally likely to be 
found anywhere on the energy shell, i.e. pE{£,t) ~ 1/T,{E) as t — oo. Inserting this into 
and applying the method of stationary phase yields 

and so the leading order term as h ^ vanishes. This implies that one has to take into account 
the next order terms. A closer analysis of shows that the important term to consider is 
the next to leading order behaviour of pE{£,t) as t ^ oo. Quite surprisingly, the two-loop 
contribution does not originate from the ergodic limit of the probability density pE{£,t) but 
from the approach to this limit. 

We have to consider pE{£,t) in more detail. It is a classical transition probability density 
and can be expressed in terms of classical trajectories. These trajectories are all time t loops 
with opening angle e. Consider one such loop as shown in Fig. Every point in the vicinity of 
its starting point is the starting point of another loop, one example being shown by the dashed 
line. To determine how angle £ and time t change with the initial point we introduce a local 
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Figure 2: a) A loop with opening angle e (full line), and the local coordinate system at its 
starting point, b) Loops with the same opening angle e whose traversal takes time t form a 
continuous family which have their starting points on a curve of constant distance d from a 
periodic orbit. 



coordinate system (see Fig. ^) and linearize the motion in the vicinity of the loop. The result 
is 

V dt = 2 cos - ds2 , V de = — 2A sin - tanh — ds2 (11) 

One finds that angle and time change only in the S2 direction, but not in the si direction. This 
is a particular property of the uniformly hyperbolic dynamics. After integrating the equations 
(|Tl]) one arrives at the following conclusion. The loops with fixed e and t form continuous 
one-parameter families. All the initial points of the loops within a family lie on a curve which 
has a constant distance (denoted by d) from a periodic orbit as shown schematically in Fig. 
The relation between the loops and the periodic orbit is given by 

cosh — sin — = cosh ^-^ (12) 
2 2 2 ^ ^ 

where to is the period of the periodic orbit. It is a remarkable property that any loop is 
uniquely related to a periodic orbit into which it can be continuously deformed through a series 
of other loops. Put another way, this implies that any self-intersection of any arbitrary classical 
trajectory is uniquely related to a periodic orbit, because a self-intersection is the initial point 
of a loop. 

We examined this property numerically. We chose a large number of long random tra- 
jectories on a Riemann surface with constant negative curvature |TI| and recorded all their 
self-intersections. For every self-intersection a point is plotted in the {e, t)-plane, where e and 
t are the opening angle and traversal time of the corresponding loop. The result is shown in 
Fig. 3. As expected, the points form continuous lines that start at the periods of the periodic 
orbits (the t- values at e = tt). One can observe a logarithmic divergence of the curves at e = 
which is implied by Eq. (12). The full line in Fig. 3 is an evaluation of Eq. (12) for the second 
family of loops, and it is found to be in perfect agreement with the numerical result. 

We continue by expressing pE{e,t) in terms of the classical trajectories. By evaluating the 
integrals over the delta-functions in (^, pE{e,t) can be written as a sum over all families of 
loops with opening angle e, which are labelled by ^ in the following. Alternatively, by using 
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7r/4 7r/2 37r/4 tt 

Figure 3: Numerical result of the search for loops with opening angle e and time t. Grey scales 
are proportional to the number of loops found in bins in the (e, t)-plane. 

the relation (12), pE{s,t) can also be expressed in terms of the periodic orbits labelled by ^o- 

. s _ _j_ ^ Tg„cosh(Ad/t,),5(t-Tg) 

PE{£,1) — Y,{E) sin|£/2|(TrM;;-2) 



S(E) ^(TrM5p-2)(TrM5p-2+4cos2 |) 

where and M^^ denote the stabihty matrices. We remark that a further use of Eq. (12) 
yields 

PE{e,t)^PE{Tr,to) . (14) 

This means that the distribution p{e, t) is identical to the return probability density p(7r, ^o) a-t 
a shifted time to^ the relation between t and to being given by Eq. (12). 

Eq. (13) is now applied to find the next to leading order behaviour of the time integral over 
Pe(£, t) as i — > oo. We assume that from a certain time Tq^e) on we can replace Pe(£, t) by its 
ergodic limit {2T^mA)~^. This time T'o(£) is chosen to have the same ^-dependence as the time 
of the families of loops (like, for example, the dashed line in Fig. 3). Thus To(£) is related to 
To(7r) by an equation identical to that between t and to (Eq- (12)). For t < To{e) we replace 
PE{s,t) by its exact form, Eq. (13). The approximation can be made asymptotically exact by 
letting To(7r) — > oo as T — > oo. We find 
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T/2 

dt pE{e, t) ~ /^/')dt ^ + ET,„<r„{.) BU^) 



+ const + 0{s^) (15) 

where here and in the following constant denotes independence of £. In the semiclassical limit 
only the asymptotic behaviour of Eq. (15) as e — is relevant and from the analog of Eq. (12) 
we find Tq{s) ~ — |log£ + const. This logarithmic divergence can be interpreted as follows. 
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For small e the two legs of a loop need a certain minimal time in order to separate enough to 
enable the loop to close. This time can be estimated by requiring that £:exp(At/2) is of order 
one, yielding the logarithmic dependence above. Substitution into (||) results in 

r^(2)/ \ 8p^T^ T-) r°°j 2e(loge+const.) 

Re j;^de' e''" e' \og{e') . (16) 

Evaluating the real part of the last integral finally yields K^'^\t) ~ — 2r^ in agreement with 
the r^-term of the GOE form factor in (^. 

In conclusion, we have shown that the off-diagonal contributions to the spectral form factor 
from two-loop orbits yield a r^-term in agreement with random matrix theory. Its origin can be 
traced to properties of loops with small opening angle e. It is expected that higher-order terms 
in the expansion of K{t) are related to multi-loop orbits, a point which is under investigation. 
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